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Abstract 

We study the asymptotic behavior of solutions to the Schrodinger equation with 
large-amplitude, highly oscillatory, random potential. In dimension d < m, where m is 
the order of the leading operator in the Schrodinger equation, we construct the heteroge- 
neous solution by using a Duhamcl expansion and prove that it converges in distribution, 
as the correlation length e goes to 0, to the solution of a stochastic differential equation, 
whose solution is represented as a sum of iterated Stratonovich integral, over the space 
C([0, +oo), <S'). The uniqueness of the limiting solution in a dense space of L 2 (ft x M. d ) 
is shown by verifying the property of conservation of mass for the Schrodinger equation. 
In dimension d > m, the solution to the Schrodinger equation is shown to converge in 
L 2 (fl x M. d ) to a deterministic Schrodinger solution in [6]. 



1 Introduction 

We consider the following Schrodinger equation in dimension d < m: 

+ ( P ( D ) ~ ^l(^)))u e (t,x) =0, t > 0, x € R d 

u e (0,x) = u (x), x <G R d , 



1.1) 



where P(D) is the pseudo-differential operator with symbol = |£| m . Taking Fourier 
transform of both sides of (1.1), we obtain 



u E (0,0 =u (0- 



1.2) 



We assume the Fourier transform of the covariance of the potential R(£) is bounded and 
continuous at 0, and the initial condition satisfies (l + |£| 2m )|^o(£)l — C uniformly in £ € M. d . 

The main objective of this paper is to construct a solution to the above equation in 
L 2 (0 x M. d ) uniformly in time on bounded intervals and to show that the solution converges 
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in distribution as e — > to the unique solution of the following stochastic partial differential 
equation (SPDE) 



i^ + P(D)u-auoW = 0, t>0,xeR d 



u(0,x) = uq(x), x £ R d , 



(1.3) 



where W denotes spatial white noise, o denotes the Stratonovich product and a is defined 

as 

a 2 := (2Tr) d R(0) = {2^) d I R{x)dx. (1.4) 

Jm. d 

To study this solution, we may take the Fourier transform of the equation: 



c u(s, x) o dW{x) 



= o{2n)- d j e~* x J e^ lX u{s,Ci)d^odW(x), 
with initial condition u(0,£) = uq(£). To look for a mild solution, we recast (1.5) as 



(1.5) 



e it m {t-s) e -i& 



u = ia{2ir)~ d 
Define formally the stochastic integral 



e^ lX u(s,^)dCids o dW(x) + u (6- ( L6 ) 



e^uis^^d^dsodWix). 



1.7) 



Uu(t,0 = {-ia){2ir)- d 
We may rewrite (1.6) as 

The mild solution to (1.3) is thus defined as u(t,x) = J-~ l {u{t, £)}, where J 7-1 denotes 
inverse Fourier transform. Suppose that d < m. The following result holds for the solution 
of (1.5). 

Theorem 1.1. Suppose d < m. The series 

ti(t,0:=2> (B) (t,0 (1-9) 

n>0 

converges in the L 2 (Q, x M. d ) sense for each t > and £ € R d and is the unique solution to 
(1.5) in the space M dense in L 2 (J7 x M. d ), which is defined in Section 5, where 



=(-ia) n (27r)- 



fc=i 



n 

k=0 



1-i 







3=1 i =1 i =1 



(1.10) 
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The following theorem shows the weak convergence of u e (t,x) to u(t,x) for any t > 
and x G R d . 

Theorem 1.2. Suppose that d < m. For any integers r > 1, mi,-- - , m r > and 
£ > " " " € K^, t\, ■ ■ ■ ,t r > we have the convergence of moments 

lira E{[u B (t 1 ,^ 1 ))]^ •••[A e (tr^ (r) )r}=E{[«(ti,e«)r i ---[«(t r ,eW)]^}. (1.11) 

e— >0+ 

The finite dimensional distribution of u(t,£) is uniquely determined by its moments of 
all orders. Moreover, the family of processes {u e (t,-),t > 0} is tight, as e — > 0+, over 
C([0,+oo);«S'(R c! )). The process {u £ (t,-),t > 0} converges in law over C([0, +oo);«S'(R d )), 
as e — > 0+, to {u(t,-),t > 0}. Also, we have that in the spatial domain, the process 
{u e (t, -),t > 0} converges in law to {u(t, -),t > 0}. 

The rest of the paper is structured as follows. Section 2 gives the formal Duhamel 
solutions for both the multi-scale Schrodinger equation and the limiting SPDE in the Fourier 
domain. Section 3 demonstrates the first order moment convergence of Duhamel solutions 
u £ . Section 4 proves that u defined as the Duhamel expansion of the limiting equation is 
the space of in L? ($1 x M d ) . Section 5 generalizes the first order moment convergence proved 
in Section 3 to arbitrary orders. Section 6 follows the approach as in [5] to show the weak 
convergence of {u e (t,£)} in C([0, +oo), S') to u(t,£) follows from tightness and convergence 
in finite dimensional distribution. 

The treatment of Schrodinger equation (1.1) in dimension d > m with potential with 
corresponding amplitude is presented in [6] . It is then shown that u £ converges in L 2 (Oxl* 1 ) 
to the solution of a homogenized equation. The asymptotic theory of solution to parabolic 
equation with large potential is presented in [1, 2]. 



2 Duhamel Expansion 

Iteratively using Duhamel's formula we obtain 



u £ (t,0 = ^ mt u (0-i^ f i e^-°)q{Qu £ (s,t-e- l Odsdt 



(2.1) 



= £4")(U), 

n=0 

where 

. n+l 



n W (U)=Hre -¥ [... [ dsy-ds n [■■■ [ TT e *(**-i-*)[€-^- 1 S^ 1 1 &r 

J J An {t) J J k=l 

77. n. \ ' J 



n$(&)#kflo(£-e -1 Z)&)- 

k=l j=l 



Here, we introduce the notation Y2j=iCj := 0) an d A„(t) := [t > si > ■■■ > s n > 0]. 
Let A n (t) := Ej=i T j — — 0]- Changing variables Sj := YH=j T i an< ^ denoting tq := 
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t — Y27=l Ti we can rewr ite (2.2) in the form 



n+l 



fin • • • dr n 



n 



,*Tfc-l 



n(t) 



fc=i 



j=i ^ i 



\ 



fe=i 

n+l 



3=1 



r+oo r+oo f f n n 

W ■ ■ / • • • / dr ■ ■ ■ dr n 5(t - Vr,-) TT g(&Kk 

Jo J J , =0 fc=1 



(2.3) 



fc=l 



Using <5(i) = J e l/3t d(3, we obtain for any rj > 

»+oo Z'+oo 



./o 

n+l 



(2.4) 



fc=i 



i=i 



Integrating out all Tj and choosing 7] = 1 we get 



it 



i B) (*,0 =H) n 



e 2 e 



/ / ••• / f[g(&K* 

J J £=i 



n 

fc=0 



1-i 



/3 



(2.5) 



fio(£-e _1 5>)- 
j'=i 



We now come to the analysis of the limiting equation. By Duhamel's principle, the 
solution to (1.5) formally satisfies the equation 



u(t,0 =e*" ,t u (0 + {-i<j){2n)- d J e ^~ s ^ 



-e* *tio(0 + (-^)(27r) - / e 

lo 



(2.6) 



e l « ia: 'u(s,^i)^idsocZW(x) 



Integrating (2.6) iteratively, we obtain formally the Duhamel expansion for (1.5) 

oo 

= E fi(n) ' ( 2 - 7 ) 

where 



n=0 



it 



(0) = e i*e m ij (^) ) and fiW = n n u {0 \ 



(2.8) 



for n = 0, 1, • • • , or more explicitly, 

=(-ia) n (27r)- nd 



A„(t) 



n+l 

n 

fc=i 



,i( s fc-l- s fc)£fc-i 



n 

i=i 



"o 



(2.9) 



fc=i fe=i 



j"=i 



Here, we define so = i, s„4-i = 0, and £o := £. At this point neither the iterative Stratonovich 
integral of u^ n ' for n = 1, 2, • • • in (2.9) nor the sum of (2.7) are well defined. We give a 
justification for these expressions in section 4. 

Using the change of variables — > £k-i ~ £k, and applying the same type of transform as 
in (2.5), we obtain 



u 



(n) 



--(-ia) n (27r)- nd e t 



/n 







k 






[-( 




'-)] 











k=l 
1 



n ^ x] Mi - e & ° n (/ir! -0'- 

j=i j=i j=i 

(2.10) 



which is the same as (1.10). 



3 Convergence of the first moments 

We shall prove the convergence of moments as stated in (1.11). For simplicity, we shall first 
consider the case when n = m = 1 and show that the limit 

lim Eu e (t,0 (3.1) 

exists. The proof of convergence for general moments is given in Section 6. 
Taking expectation of both sides of (2.1), we obtain 

Eu e (i,0 = E E ^ 2B) ( t '0- (3-2) 

n>0 

This is because expectation of product of odd number of Gaussian random variables is 0. 
The expectation of product of even number of Gaussian random variables is given as a 
sum of products of the expectation of pairs of variables, where the summation runs over all 
possible pairs. The contribution of products of potentials can thus be represented by 

in 

E{n<z(£ fe )}=E n (3.3) 

fc=l 7T (e/)67T 

where (ef) denotes pair of indices, ir denotes a pairing of the 2n indices, and the summation 
is over all possible pairings. 
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Adding up all the delta functions gives $Zfc=i 6c = 0. We can therefore write 



( 2n 










x n 


['-'( 


e-^E^ 


'-): 








3=1 







2/i 



fc=i 



<-i)" £ ^ (e)E^ e ^^ /•••/ 



f 2n 










n 


[,..( 


e-^E^ 


'-): 


i 






i=i 







[] i?(e e )5(e e + ^)de e ^/- 

(e/)er 

(3.4) 

The above summation extends over all possible pairings 7r made over vertices {1, • • • , 2n}. 
By changing variables := we obtain 



En( 2 ")(t,O=(-l)Vn (OE J ( e ^ (i/3 / "' / 



C 2n 




k 


m 


x n 


[-( 


e-Eo 








3=1 





|8 



(e/)er 



In what follows, we show that 



sup|E4 2n )(t,e)| <a n , 



where ^ n a n < +oo. Hence, 



lim EuJt,£) = V lim Eui 2n \t,0. 



£->0 + Z — ' £— >0-| 

n=0 



(3.5) 
(3.6) 

(3.7) 



Denote by C(n) the set of all left vertices of a given pairing it. Let 

fe 



Ai 



3=1 



Define 



= [l-^-/3)]- 1 . 
Using contour integration method, we are able to show that 

(p) _ / e iAt tP- l e- 1 t > 0. 
&A ~ \ f < 

The following estimate is then derived: 



(3.« 



(3.9) 



(3.10) 



Lemma 3.1. Suppose that p > 0. There exists a constant C > such that for an arbitrary 
n > 1 we have 



e w\l[[l-i(A k -f3)]\ d(3 



.k=l 



< 



C n t np-l e -t 

[{n-lW ' 



(3.11) 



Readers are referred to [5] for the proof of this lemma. 
With such notation, we may rewrite (3.5) in the form 



Eu^(t,0 =(-l)Vfio(OX; /••• / II *e Ak *F{t,^) 

T k(££{n),k^2n 

x n fl( e e e )^ e +e/)de e de/, 

(e/)67r 



(3.12) 



where 







k 






Kc-P))- 2 n 


h( 




"-«): 


d/3 


fce£(7r) 











d/3 (3.13) 



an d rifc^£(7r) k^2n *&A k denotes the convolution of all e^ fc 's, where is the right endpoint in 
the giving pairing ir, except k = 2n. 



Lemma 3.2. Suppose p € (—, 1). Then 

[ | -p+ \£-u\™ + i\- p d£, < +oo. 

•Jm. d 



sup 



(3.14) 



Proof. We may first shift £ to get rid of w and then change variables to polar coordinates 

r r+oo 

/ \-p + \z-u\ m +i\-f> <n d / i -/s+r+^r^-^iei, (3.15) 

where is the area of the unit sphere in M. d . Let Q = £ m . The integral on the right hand 
side of (3.15) may be rewritten as 



r+oo 

/ \-P + Q + i\~PQ d / xa - 1 dQ. 
Jo 



(3.16) 



Without loss of generality, we assume (3 > and the integral above can then be written as 
the sum of three integrals /, II, III according to whether £ belongs to (0, /3/2), (/J/2, 2/3), 
or (2/3, +oo). We have 



'o 

The second integral can also be estimated as follows 



II<(2^~ 1 I \Q-p + i\~PdQ. 



(3.17) 



(3.18) 
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If (3 < 1 we have 

If /3 > 1, we estimate 



IK 3 x 2i" 2 . 



// < (2/3)^" 1 y |Q - /3\-PdQ < [1 + (-) 1 ~12^- 
The third integral is estimated 



/JJ<^ + l)-5^<^. 



This concludes the proof of the lemma. 
Let 



(3.19) 
(3.20) 

(3.21) 
□ 



Fx 



:= / 

Jr 



[i - ice 1 - /3)] 2 













• n 


[„.( 


e-E& 


'-): 


1 


fce£(7r) 




3=1 







(3.22) 



and 



F 2 (U;7r) := / e^d/3< 







fc 






n 


[-( 


e-E<* 


■-): 


-(i-p) | 


fc6£(7r) 




3=1 







Since 



F 1 (t^; 7 r) = e ^ ) *...*eg +2(0 (t) 



(3.23) 



(3.24) 



for some Bj(£), we get Fi(i,£;7r) = for i < 0. Likewise, F2(t,£;ir) = for t < 0. We can 
therefore write 



F(i,£;7r)= I F 1 (t-s,&7r)F 2 (s,Z;7r)ds>0. 
Jo 

Observe that by Lemma 3.1 

(jn^n(l— p) — 1 



F 2 (t,£;iv) < 



and 



n *^ k {t)< 

for i > 0. On the other hand, we have obviously 

dp 



C n t n-1 



(n-1)! 



(3.25) 

(3.26) 
(3.27) 



|*l(t,fc7r)|<C7/ - 



+ IC m -/5| 2 



f 




k 






n 


-( 




'-) 


1 


1 fcg£(7r) 




3=1 







(3.28) 
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for some constant C > 0. As a result, we obtain 

"* (t - sj^ds 



|E4 2n) (^e)i<ce*K(oiE 



o (n-1)! 



(W)6tt 



< 



CV|fi (Q| 

[( n _l)!]2-p 















-( 




"-) 




1 fce£(7r) 











(e/)67r 

(3.29) 

Here, we just need half of the order of decay of the initial condition «o(£) as assumed to 
have 

< C 1 

(1 + |£™ - /3| 2 )V2 - (1 + £2m)l/2 (1 + |£m _ 0(2)1/2 ~ (1 + 02)1/2 " V • ; 

The right hand side of the above inequality is then used for estimating the integration in (3 
by Lemma 5.2 in [6]: 



+oo 



l + log+|g| 

(1 + /3 2 )V2(1 + |£m _ 0(2)1/2 - ~ (1 + ^2m)l/2 ' 



< c 



(3.31) 



where log + |£| := max(0, log |£|). 

Using Lemma 3.2 we conclude therefore that 

|E4 2n) (U)l<(2n-l)!! 



C n t n{2-p) e t l + log + |f| 
( n !)2-p (l + £2m)l/2 



(3.32) 



and (3.6) follows by letting d2 n = (2n - l)!! c " ( ^,_p )et . We obtain 



(n\) 2 -P 
+oo 



(3.33) 



n=0 



where 



^)(t,e) H-£(0))Vfi (£)E / / 

2n 

n 



fc=0 



1 - i 



5=1 



/3 



-i 



(3.34) 



n <y(6+e/Kede/. 

(e/)67r 



In what follows, we show that 



t (2n) (t,£) =En( 2? %£). 



(3.35) 
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Expectation of multiple Stratonovich integrals is denned: 

2n 

E iYl dW ( x j)} = E II 6 ( Xe ~ x f) dx edXf. 
3=1 i" (e/)G7r 



(3.36) 



Upon integrating in all variables x, the first moment of u^ 2n ^ can therefore be written as 



2« 



Eu^ =(-ia) 2n (27r)- 2nd e t V / • • • / / / TT ^ 



( 2n 




k 






M 




■-.)] 











(3.37) 



mo n ^+^' 

(e/)€7r 



which gives exactly (3.35). 

4 I? convergence of the SPDE solution 

We now come to the series (1.9), which we claim is the solution to (1.6). The readers are 
referred to [1] for more complete description of the relationship between Stratonovich and 
ltd integrals and for additional details on the theory presented in this section. 

We first prove that u(t, •) in (1.9) as a series is well-defined in the space of L 2 (fi x M d ). 
Denote = I n (f n ), where l n denotes the n-th order iterated Stratonovich integral, and 
f n is a n-parameter function, i.e. 



-E-n{fn) 



f n (xi, x n )dW(xi) . . . dW{x n ). 



(4.1) 



By the definition of the L 2 norm of multiple Stratonovich integral, we have 

E\u^\ 2 = E{l 2n (f n ®f n )} 

= E I f n (x 1 ,...,x n )f n (y 1 ,...,y n )dW(x 1 )...dW(x n )dW(y 1 )...dW(y n ). 

(4.2) 

Taking the same approach as in the calculation of the first moment in (3.37), we obtain 



E\u^(t,0\ 2 = (a n (27rr d e t ) 2 J2 [ [ [ e^^d^dhMi ~ £ S 3 



T2n 



s.(f+Eff)x n 

j=l [fc= 

II ^e+£/K (1) ^ (2) . 

(e/)G7r 



l-i 



f-E< 



(i) 



i=i 



n 

k=0 



l+i 



(2) 



i=i 



+ /3 2 



(4.3) 
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We denote d£® = d^p ■ ■ ■ d£$ ,1 = 1,2. £ e and are paired arguments in the graph 
comprised of 2n arguments in total with their index e and /. 

The idea that we use to estimate the first moment applies here too. However, note that 
we have {k = n} G £(7r) for the crossing graphs defined in [2], and one of the terms [1 — 
ii^-P)]- 1 in the function F x defined in (3.22) now becomes [1 - YJj=i tf* T 
Therefore, we have to make the following adjustments in the proof. Using the smoothness 
condition for the initial condition, we have 



i-i(K-£Li€ 



(1) i™-/v 



< 



(i + ic-E"=ie 7 (1) i 2m ) 1/2 ( 1 + /3 i) 1/2 ' 



(4.4) 



We use the term [1+|£— Z)j=i Cj | 2m ] ^ 2 f° r tne integration in £„, and the term [1+/3 2 ] x l 2 
together with [1 — z(£ m + /3i)] _1 indexed by {k = 0} in the first product in (4.3) 



1 + log+ j|j 

(1 + /32)1/2|! _ i^ra + ^)| -= " (1 + £2m)l/2 " 



< c 



(4.5) 



We also use the smoothness condition for the initial condition and obtain another (1 + 
log + |£l)/(l +£ 2m ) 1 / 2 from the integration in fa- Finally, we obtain the estimate 



E\u^ n \t,0\ 2 < (2ra-l)H 



C n t n(2-p) e t l + log + |£| 



(n!) 



(l + £2m)l/2 



(4.6) 



for any p € (^,1). This implies that the iterated Stratonovich integral u^ n ' is indeed 
well-defined. Integrating in £ and summing the above bound over n gives 



E I \u(t,0\ 2 d^ = J2i E [ \u {n \t,t)\ 2 d0 1/2 

n>0 ^ 



< OO 



(4.7) 



and the L 2 (Q x M, d ) convergence of (1.9) follows. In fact, by first multiplying £ m with (4.6) 
and performing the integration and summation we can futher obtain 



E 



J CWt,0\ 2 dZ<oc. 



(4.8) 



5 Uniqueness of the SPDE solution 

Let us now provide the rigorous definition for the operator %. Suppose f(t,£,) is a sum of 
iterated Stratonovich integrals 

/(t,0 = Z) £,-))■ (5.1) 

n>0 

We define 

n>l 
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where 



(Hf)n+i(t, , x n) ■= (-i<r)(27r)" 



/n(s,£o,zi, • • • ,x n )d£ods. 

(5.3) 

It is checked that under this definition Tiu = X^n>i ^ = u — u^°\ The Duhamel solution 
(1.9) is therefore a solution to the equation (1.6). 

We can also define 

Jf(t,0 =^X„+i((J7)n+l(i, £,0), (5.4) 



n>0 



where 



(J7)n+l(*,£,z,xi,--- ,z n ) := (-ia)(2n)- d e-^ x J e*> x f n ( a , faxi, ■ ■ ■ ,x n )d£ . (5.5) 
For the Duhamel solution (1.9), we have 

n+l 

E|X„ +1 (Un) n+1 (t,0)| 2 = K(2vr)-V) 2 ^ / / f e^-^d^dfoMt ~ £ ff) 



n+l [ n+l 

xs.«+E«f) n 

j=l ^ fe=l 

II *(&+e/)de w de (3) 

(e/)G7r 



1 - t 



(i) 



+ /?i 



n+l 

n 

it=i 



i + t 



(2) 



i=i 



+ /3 2 



Although it looks a little different from (4.3), it can be estimated in the same way as 



(5.6) 



(E|«( n )| 2 )V2 < (2n- 1)!! 



C n t n{2-p) e t 1 + l og+ |£| 
( n !)2-p ( 1 + £2m)l/2> 



(5.7) 



which uppon summation in n implies that E|n(t,^)| 2 is uniformly bounded for all £ € 

f.' 

the equation 



Note that Uu{t,£) = (-icr)(2vr)- d J * e^*"^™ {Ju){s, £)ds. u(t,£) is therefore a solution to 



(^+rHi,0 = ^(2vr)- d Jn(t,e). 



(5.8) 



Now we prove that this equation preserves mass. Multiplying this equation by u(t,£), and 
integrating in £ and over the probability space gives 



~E / \u\ 2 d£ + E I C\u\ 2 d£ 



(5.9) 
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The right hand side of this equation can be written out explicitly as 

E J (Ju)~udi = J ]TE (e-* x J e^ x f n (t, 6)^1) Zm 0) } ^ 



/E / / ei?lX /n(^ei)^l/m(t,C)E II ^e-^)dx) d£ 

EE / e ^ 1X fn(t, HlHl j e-^f m (t,i)di J] J(a; e 
n,m tt J (ef)Ew 

EE/ fn(t, x )fm(t,x) Y[ 5(x e -Xf)dx, 



n,m ir 



(e/)67r 

(5.10) 

which is real-valued because of the symmetricity of this summation. Extracting the imagi- 
nary part from both sides of (5.9) gives 

d 



Ot 



E j \u(t,0\ 2 d£ = 0. (5.11) 



Finally, we define the space M in which the equation (5.8) admits a unique solution. In 
light of the equation (5.9), M consists of sum of iterated Stratonovich integrals f(t,£) = 

En>0 2 "(/«(*'^' ')) such that 



1. f(t,0eL 2 (fix 

2. Jf € L 2 (n), 

3. |#/(t,£)€L 2 (Ox 

As a reminder, defining the sum of iterated Stratonovich integral / as in (5.1), we have 
that 

n>0 m>0 

where 

^ ( m _|_ 2k)\ f 

9m(t,£,x) = 22 | fc , 2 fc / fm+2k(t,£,x m ,yf 2 )dy. (5.13) 
Here, y <8>y = (y,y), and 

I'm \9m ) '■= / 9m{t,£,xi,- ■ ■ ,x m )dW(xi) ■ ■ ■ dW(x m ) (5-14) 

JR md 

denotes the iterated Ito integral. The L 2 norm of / can then be computed using the 
orthogonality of Weiner Chaos expansion as 



2 x 1 



l/IU»<n> = E™ ! / E^^LlWI^^^yf^y dx <oo, 

\m>0 17 \fc>0 ' • J) 

(5.15) 



13 



The readers are referred to [1] for more details. 

It is easy to verify that the space defined above is dense in L 2 (J7 x Denote the space 
consisting of all functions that satisfy condition (1) and (2) by M. In fact, any function 
/(£) £ L 2 (Q x M. d ) can be written as their Wiener Chaos expansion 

/(£) = X;JmGU£,Xm))- (5-16) 

m>0 

Each g m can be approximated by a function which vanishes in a set of measure at most 
k^ 1 in the vicinity of measure set of diagonals given by the support of the distributions 
5(x e — Xf). By the change of change of coordinates in (5.13), we have gm = fm 1 so that 
the ltd and Stratonovich integrals agree. Define 

/W(e) = ^X m (/«). (5.17) 

m>0 

Using formula (5.15), we may verify that 

\\Jf ik) \\mn)<^ (5-18) 

and 

Urn \\f {k) (O-f(O\\mnxR*)=0- (5-19) 

We have shown that M is dense in L 2 (f2 x R rf ). Since M is dense in M, it is also dense in 
L 2 (n x R d ). 

6 General moment convergence 

We now turn to the general case of (1.11). It suffices to consider the limit 

limE{u e (t 1 ,^)---u e (tr,& ) )}- (6.1) 
Using (2.1), we obtain that the above expression before passing to the limit equals 

£> £ (n), (6.2) 

where 

X e {n) :=E{u^(t 1 ,^)---4 n ^(t r ,e ) )}. (6.3) 

The summation in (6.2) extends over all non-negative integer multi-indices n = (nx, • • • , n r ). 
We adopt the customary notation \n\ = ^[=i n i- I 11 fact, only the terms for which |n| are 
even do not vanish. 
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Each term appearing on the right-hand side can be written as 

r „ r r „ „ r 

l £ (n) :=H)H ex p(V^) V / eMiyZmU^i / / TT^ 



"I 



i=l fe=l (e/)G7r 



(6.4) 



x < 



r n\ 










nn 


[-( 




"■')] 


r 


1=1 k=0 

V 




3=1 







We denote := • • • d£n t . The summation extends over all possible pairings. We can 
repeat the argument made in the previous section and obtain bounds of the form 



|2e(n)| < (H-l)!! 



{Ct 2 -P)^e TT ) 



(6.5) 



assuming < t%, . . . ,t r < T, where the constant C does not depend on either e G (0, 1], 
|n| or r. There are (^^T^ -1 ) non-negative integer valued multi-indices satisfying equation 
In I = N. We can therefore estimate 



E u< 

n\=N 



< CAT, 



(6.6) 



where 



cat := 



r-l ) K ' (f!)2-P 



(6.7) 



and the series ^at=i c iv is clearly summable. We are therefore allowed to pass (6.2) to the 
limit e — > 0+. As a result, we conclude that 



lim E{u e (t 1 ,H^)---u e (t r ,^)} = J2 lim X e (n) = Vx(n), 



(6.8) 



where 



Z(n) :=(-ia(2nr d ) lnl Yl [ exp{iJ2 mfldfr [■■■ f f[d£® 

n jRlnd \ 1=1 1=1 J J 1=1 



"I 



1=1 k=l (e/)67r 

f r ni 



>< nn 

^=1 fc=0 



1 -t 



3=1 



(0 



(6.9) 
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7 Weak convergence 



To show weak convergence, we need to first demonstrate the convergence of finite dimen- 
sional distributions, i.e., the convergence of distributions of (u e (ti,£W), ■ ■ ■ ,u E (t r ,£^)) for 
an arbitrary r > 1, t\, ■ ■ ■ ,t r > 0, • • • , £( r ) G M d . Since we have already proved the 
moment convergence in the previous sections, it suffices to verify the determinacy of the 
distributions by their moments. For simplicity we consider only the case r = 1. It is easy 
to generalize to the case for arbitrary r. Using estimates (6.6) we obtain that 

N=0 V 7 V 2 ■> 

Using Stirling's formula we can easily obtain that 

N + r ; 1 )<C(l+ r -^f( 1 + ^-y-K (7.2) 
r — 1 J TV r — 1 

Plugging this into equation (7.1) therefore gives 



|E[«(t,e)n<C2 r -V*^ (l + -— (AT-l)!! 

at— n \ / 



N=0 v 7 v 2 

±2S / iV V"* ,„(Ct 2 -")fe rt 



N=r V r / 



(7.3) 



The first term can be easily estimated by Cr r while the second by a constant C independent 
of r. Therefore we have 



-oc 



r=l L v /J ri=l 

and the uniqueness follows from Carleman's condition. 
It remains to prove the tightness of {u e (t, £)} over 

C([0,+oo),«S'(R d )). By [4], it suffices 
to prove that {uf(t) := (u £ (t, ■),(/)) L 2n^d\,t > 0} for an arbitrary (ft € S(M. d ), which, by 
Kolmogorov's theorem, follows from: 



Proposition 7.1. For any T > and (ft G <S(R ) there exists a constant C > swc/i that 

E\uf(t) -uf(s)\ 2 < C(t-s) 2 ,Ve € (0,1], s,t G [0,T]. (7.5) 
Proof. From equation (1.2), we obtain 

4(t) - 4{s) = if T u*(r)dr + ^(r)dr, (7.6) 
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and 



( T ) :=J>^(t) 



n>0 



d(n + l) 



(7.7) 



x < 



n 

k=0 



1-i 



i=o 



■/3 



-i 



tio(£-e _1 Z)W(0- 

5=0 



As in section 4, we conclude that for all r € [0, T], we have E|i£ m iie (t)| 2 < C and E|ue (r)| 2 < 
C for some constant independent of e € (0,1]. Estimate (7.5) is a consequence of the 
Cauchy-Schwarz inequality. □ 

Using (2.10) for any <fi € <S(IR d ), we can write 



e T dr / e 



n 

fc=0 



1-i 



fc=i 

m 



+ P 



-1 



"0 



5=1 



5=1 



5=1 



■t y e r dr / ^ 

n 

n 



k=X 



1-i 



*-E& +^ f ^o(e-E^)II e "^ o II d ^)- 

5=1 / J J 5=1 5=1 5=1 

(7.8) 

Applying the same technique as in the proof for the L 2 (]R rf x O) setting, we check that for 
any T > 0, we have 

E|(nW(V),0) - (nW( S ,-),0)| 2 < (CT) 2 (i- S ) 2 (n!) (1 ~ p) ,V™ > l,s,t e [0,T] (7.9) 

for some constant C > independent of n. This in turn implies that 

E\(u(t, ■),(/)) - (u(s, •),</>) | 2 < C(t-s) 2 ,Vn > 1 (7.10) 

on any compact set. By the Kolmogorov-Chentsov Theorem, we have that (u(t, ■),([)) is 
continuous almost surely. 

Weak convergence of {u e (t,£)} follows from the convergence of finite dimensional distri- 
bution and tightness. Returning to the spatial space, by the Plancherel theorem, we have 
(u(t, ■),(/>) = (u(t, •),(/>)■ Hence, the process {u e (t, x)} converges in law over C([0, +oo); S' 
to {u(t, x)}. 
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